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We propose a continuum representation of the Dynamical Mean Field Theory, in which we were 
able to derive an exact overlap between the Dynamical Mean Field Theory and band structure 
methods, such as the Density Functional Theory. The implementation of this exact double-counting 
shows improved agreement between theory and experiment in several correlated solids, such as the 
transition metal oxides and lanthanides. Previously introduced nominal double-counting is in much 
better agreement with the exact double-counting than most widely used fully localized limit formula. 
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Understanding the electronic structure of materials 
with strong electronic correlations remains one of the 
great challenges of modern condensed matter physics. 
The first step towards calculating the electronic struc¬ 
ture of solids has been achieved by obtaining the single¬ 
particle band dispersion E{ k) within the density func¬ 
tional theory (DFT) in the local density approximation 
(lda) n, which takes into account correlation effects 
only in a limited extent. 

To account for the many-body correlation effects be¬ 
yond LDA, more sophisticated methods have been devel¬ 
oped. Among them, one of the most successful schemes is 
the dynamical mean-field theory (DMFT) [2]. It replaces 
the problem of describing correlation effects in a periodic 
lattice by a strongly interacting impurity coupled to a 
self-consistent bath |2j. This method was first developed 
to solve the Hubbard model, but it was soon realized [4] 
that it can also be combined with the LDA method, to 
give more material specific predictions of correlation ef¬ 
fects in solids. The LDA+DMFT method achieved great 
success in the past two decades, as it was successfully 
applied to numerous correlated solids [5]. The combina¬ 
tion of the two methods, nevertheless lead to a problem 
of somewhat ambiguous way of subtracting the part of 
correlations which are accounted for by both methods. 

The so-called double-counting (DC) term was usu¬ 
ally approximated by the formula first developed in the 
context of LDA+U, and was evaluated by taking the 
atomic limit for the Hubbard interaction term 0 m. 
Many other similar schemes were proposed recently 0- 
m, but rigorous derivation of this double-counted inter¬ 
action within DMFT and LDA is missing to date. Here 
we propose a new method of calculating the overlap be¬ 
tween DMFT and a band-structure method (either DFT 
or GW), and we explicitly evaluate this DC functional 
within LDA+DMFT applied to well studied transition 
metal oxides such as SrV 03 , LaV 03 , and most studied 
lanthanide metal, the elemental Cerium. 

To compare different approximations in the same lan¬ 
guage, it is useful to cast them into the form of the 
Luttinger Ward functional Banana, which is a func¬ 


tional of the electron Green’s function G, and takes the 
form r[G] = -Tr((Go 1 -G- 1 )G) + Trlog(-G) + $y c [G]. 
The first part is the material dependent part, in which 
Gq 1 (rr / ; u) = [u> + p + V 2 - V ext (r)]S(r - r'), and 
the second two terms are universal functionals of the 
Green’s function G(rr, v't') and the Coulomb interac¬ 
tion V c {r — r'). In the exact theory, $y c [G] contains all 
skeleton Feynman diagram, constructed by G and V c Ca¬ 
in the language of Luttinger Ward functional, different 
approximations can then be looked at as different ap¬ 
proximations to the interacting part <I>y e [G]. 

The Density Functional Theory can be derived by ap¬ 
proximating the exact functional $y c [G] by E H \p(r)\ + 
E X c[p{ r )], where Eh and E xc are the Hartree and 
the exchange-correlation functionals, respectively. The 
stationarity condition gives the DFT equations, i.e., 
G -1 - Gq 1 = (V H [p\ + V xc [p])5{v - r ')8(t - r'), be¬ 
cause 6E xc [p]/SG = <5(r - r')5(r - +) 5E xc [p\/Sp = 
S(r — r')S(T — t') V xc [p). Note that in this language, exact 
DFT appears as an approximation to the exact Green’s 
function, where the exact self-energy is approximated by 
a static and local potential. Note also that the static 
approximation is a consequence of truncating the vari¬ 
able of interest, namely replacing full G(r, r, r', t') by its 
diagonal components p( r) = 5(r — r')<5(r — r')G(rr, r'r'). 

In the Luttinger-Ward functional language, the DMFT 
appears as an approximation where the Green’s function 
in the <f> functional is replaced by its local counterpart 
G -+ Gi oca i , and the Coulomb repulsion V c by screened 
interaction V c -A 17, namely $ DMFT = d>(y [G; oca ;]. 0 
Note that the DMFT functional has exactly the same 
form as the exact functional <f>y c [G], because all the 
skeleton Feynman diagrams constructed by Gi oca i and 
U are summed up by DMFT, while in DFT the func¬ 
tional E xc [p] is unknown, and further approximation is 
necessary. The truncation of the variable of interest from 
G to Gi oca i leads in DMFT to self-energy, which is also 
local in space, but it keeps its dynamic nature. Other 
approximations such as Hartree-Fock or GW can be sim¬ 
ilarly derived by replacing $y c [G] by some limited set 
of Feynman diagrams, i.e., truncation in space of Feyn- 
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man diagrams, rather than truncation of the variable of 
interest. 

There is some kind of disconnect between the DMFT 
functional $ FMFT [Gi oca i ], and the LDA functional 
E xc [p{ r)], mostly because the auxiliary systems for the 
two methods are very different. The auxiliary system for 
LDA approximation is the uniform electron gas problem 
defined for continuum, in the absence of complexity of the 
solid. On the other hand, DMFT is usually associated 
with the lattice model like Hubbard model, where map¬ 
ping to the local problem reduces to the Anderson impu¬ 
rity model, which does not have a well-defined continuum 
representation. The double-counting problem occurs be¬ 
cause it is not clear what is the overlap between the two 
methods, i.e., what physical processes are accounted for 
in one and what in the other method. 

It is useful to represent the DMFT method in the 
continuum representation. Such representation is not 
unique, but physical intuition can guide the map¬ 
ping. Here we propose to look at the DMFT prob¬ 
lem as the approximation, which solves exactly the 
problem defined by some auxiliary Green’s function 
Giocai = PG and Coulomb repulsion replaced by Yukawa 
short-range interaction VP = e A-pr 1 . We have in 

mind some projector P, which is very local, and trun¬ 
cates the Green’s function to a region mostly concen¬ 
trated inside the muffin-tin sphere. It can for ex¬ 
ample be defined by a set of quasi-atomic orbitals 

Giocaiir,^) = ( r l$i) ( $ il G\(/) L ') W>z/I r ') wher e 

(r|$i) = ui{t)Yl(y) are spheric harmonics times lo¬ 
calized radial wave function. Note that this trunca¬ 
tion of the Green’s function G(r,r') to its local coun¬ 
terpart parallels the truncation of the Green’s function 
to its diagonal component in theories that choose den¬ 
sity as the essential variable, i.e., p(r) = G(rr, y't’)5{jc — 
r')S(r — t'). The screening A in Yukawa interaction Vp 
has to be large enough such that the interaction be¬ 
tween electrons on neighboring sites is negligible. The 
DMFT can then give an exact Luttinger-Ward func¬ 
tional | I > yA[PG], i.e., containing all local Feynman dia¬ 
grams constructed by PG and V*, defined in the contin¬ 
uum mi. The stationarity condition for the Luttinger- 
Ward functional gives the DMFT equations G _1 —Gq 1 = 
P (£$Vc* [ Gl OC al ] / SGl oca l ). 

The precise determination of the screening A is beyond 
the scope of this paper. However, we notice that once the 
Coulomb interaction U in DMFT is known, the screening 
length A is uniquely determined by U through the ma¬ 
trix elements of the Yukawa interaction in DMFT basis. 
Notice that Hund’s coupling J is not a free parameter 
in this parametrization, but is uniquely determined by A 
through Yukawa form of the Coulomb interaction. [21] 

After the mapping of the DMFT method to the 
continuous (r,r / ) Hilbert space, where DFT exchange- 
correlation is defined, it is easy to see what is the 


overlap between the two methods. The Hartree term 
is accounted for exactly in the LDA method, and has 
the form Ey c [p\ = \ f drdr'p(r)p(r')V c (r — r'), while 
in DMFT it takes the following form e H ’ DMFT = 
\ f drdr'(Pp(r))(Pp(r'))Vp(r — r'), which can also be 
written as e h,dmft = Ey X [Pp \, where Pp = S( r — 

r')6(T-T , )G; oca ;(rr,rV) = <5(r-r , )<5(r-r , )PG(rr, rV), 
and Ey \p] is the exact Hartree functional defined above. 
The Hartree contribution to the DC within LDA+DMFT 
(or any other band structure method which includes ex¬ 
act Hartree term) is thus Ey X [Pp] [T5j. This DC term 
thus corresponds to truncating Green’s function G and 
Coulomb interaction V c by their local counterparts, i.e., 
G —)• PG and V c -+ Vp. 

For approximations, which truncate in the space 
of Feynman diagrams (such as Hartree-Fock or GW 
method), one can obtain the DMFT double-counting by 
applying both the truncation in space of Feynman dia¬ 
grams as well as the DMFT truncation in variables of 
interest. For the case of GW method, one can check di¬ 
agram by diagram that the corresponding DMFT Feyn¬ 
man diagram is obtained by replacing G by PG and V c 
by VP in each diagram, just like it was done above for the 
Hartree term. More precisely, the GW functional can be 
written as <l> { p w [G] = E F - §Trlog(l - V C G * G), where 
G*G = P is the convolution of two Green’s functions (po¬ 
larization function). The GW+DMFT double-counting 
is thus E h < dmft - ±Trlog(l - Vp{PG) * (PG)), which 
can be shortly written as ®PPP [PG]. 

* c 

In the case of DFT+DMFT, the expansion in terms 
of Feynman diagrams is not possible, however, to iden¬ 
tify the overlap between the two methods, this is not 
essential. Clearly, the double-counting in DFT+DMFT 
is obtained by the same procedure of replacing G by PG 
and V c by Vp in the DFT functional. Since the DFT 
also truncates the Green’s function to its diagonal com¬ 
ponents only (p = S(t — t')S( r — r')G) the DC is a func¬ 
tional of the local charge only pi oca i = Pp. DC thus takes 
the form 

^ DFT+DMFT = ^ { p p] + E*f[Pp\. 

In LDA method, the exchange-correlation functional is 
obtained from the energy of the uniform electron gas. To 
obtain the LDA+DMFT double-counting, one thus needs 
to solve the problem of the electron gas with the density 
that contains only ’’local” charge Pp, which interacts by 
the screened Yukawa interaction Vp. [23] 

Including the exact double-counting, the LDA+DMFT 
$ functional is thus 

^LDA+DMFT^ = + E ^ C [p] + $yx[PG} - 

-E F x [Pp\- E™[Ppl (1) 

where [PG\ is the DMFT functional which contains 
all Feynman diagrams constructed from PG and V) A . 
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This is the central equation of this paper, as it defines 
the LDA+DMFT approximation including the exact DC. 
The saddle point equations give the LDA+DMFT set of 
equations 


Ce-a 

n f 

V dc /U 

exact 

0.997 

0.424 

nominal 

1.002 

0.500 

FLL 

1.035 

0.533 


G- 1 -Gn 1 = P 


M'Fy* [ Glocal] 


SG, 


local 


(SE* xc \p\ „6E» xc [p local 

1 - _ JJ 


^ s P 


S Plocal 



TABLE I: LDA+DMFT valence and DC potential for a-Ce 
at T = 200 K. The local Coulomb repulsion in Ce is U = 6 eV. 


where 


where we used E HXC [p] = E H [p\ + E xc [p\ and PG = 
Glocal • 

The only difference between functional Eq. [TJ and the 
usual LDA+DMFT implementation, is the presence of 
E^ xc ■ This is the semi-local exchange and LDA correla- 
tion functional of the electron gas interacting by Yukawa 
interaction. The semi-local exchange-density £?+ [p] (de- 
fined by E x [p\ = J drp(r)e x [p(r)]), can be computed an¬ 
alytically, and takes the following form 

e v c M = 

' S 


where 


log(l + ai) = 


A(cro + OqA) 


1 + a^A 2 + a 3 A 4 + aqA 6 

log(1 + “ 2 ) = 1 + & v + & v 

, n , N A 3 ( 7 o+ 7 iA) 

108(1+ ° 3)= 1 + 72 A» 

log(l + 04) = A 4 (Jq + ^iA 2 ) 


( 4 ) 

( 5 ) 

( 6 ) 
( 7 ) 


The best fit gives the following coefficients: 

= [1.2238912,7.3648662,9.6044695, 

—0.7501634,0.0207808] x 10 " 1 
Pi = [5.839362,11.969474,10.156124,1.594125] x 1CT 2 
7 i = [8.27519,5.57133,17.25079] x 10 -3 
Si = [5.29134419,0.0449628225] x 1CT 4 ( 8 ) 




1 4 arctan(2x) (12x 2 + 1) log(l + 4x 2 ) 


3x 


+ 


24x 4 


c = If A) 


1/3 


2 V 47r 2 ) 5 


= 

V 4 n P ) 


1/3 


. o»^ = (¥) 1/3 +- 


The exchange potential V x = j^E x [p\ is then Vy X = 


4 x 1 1 

3 c- T/A I o ^ . 


3 VA 


3 r B dx' 


The correlation part requires solution of the homoge¬ 
neous electron gas problem interacting with Yukawa re¬ 
pulsion, which was solved by QMC [751 - 178] . Here we 
want to have an analytic expression for correlation en¬ 
ergy at arbitrary A and r s . It is well established that 
G 0 W 0 gives quite accurate correlation energy of the 
electron gas [19] [2D], especially when computed from 
the Luttinger-Ward functional F[G]. We thus repeated 
GqWq calculation for the electron gas, but here we use 
Yukawa interaction. We evaluate the total energy using 
Luttinger-Ward functional of GW to achieve high accu¬ 
racy. We then fit the correlation energy in the range 
of physically most relevant r s £ [ 0 , 10 ] and screenings 
A £ [0,3] (A is measured in Bohr radius inverse) with the 
following functional form: 


£ 


c 

w 


'A=0 


£ 4 „=1 


a n r" 


( 3 ) 


Finally, the correlation potential is V£ = ’ 

where A(r s , A) = 1 + a n r s and C(r s ,A) = 3[1 + 
£n=i a n r s ] 2 /£l=i n a n r s ■ We take the unscreened 
correlation energy density £^ =0 (and unscreened po¬ 
tential) from the standard parametrization of quantum 
Monte Carlo results, hence the G 0 W 0 calculation is only 
used for renormalization of correlations by screening with 
Yukawa form. 

In the following we present results for some of the 
most often studied correlated solids, namely, elemen¬ 
tal Cerium, SrV 03 and LaVOa. We used three differ¬ 
ent forms of DC functional: i) ’’exact”, which we intro¬ 
duced above, ii) ”FLL” stands for fully localized limit 
form introduced in Ref. El which has the simple form 
Vdc = U(n— 1/2) — J/2(n— 1), and n stands for the cor¬ 
related occupancy, c) and the ’’nominal” DC, introduced 
in Ref. [9, 10l The ’’nominal” Vd c takes the same form as 
”FLL” formula, but n in the formula is replaced by the 
closest integer value (n° = [n]), and hence n° corresponds 
to so-called nominal valence. We use LDA+DMFT im¬ 
plementation of Ref. H 

The physical properties of correlated materials are 
very sensitive to the value of the local occupancy n/, 
and nf is sensitive to the value of DC. In table [T] we 
show results for elemental Cerium in the a phase. All 
three DC functionals give very similar correlated occu¬ 
pancies rif, and all are very close to nominal valence 
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SrV03 

n t2g+eg 

n t2g 

n eg 

v dr 9 / u 

V dr./ U 

exact 

2.223 

1.507 

0.716 

1.384 

1.406 

nominal 

2.251 

1.541 

0.710 

1.443 

1.444 

FLL 

2.529 

1.699 

0.830 

1.943 

1.943 


LaV03(t2g-only) 

n t2g 

v£ g /u 

vi 9 /u 

exact 

2.014 

1.195 

1.193 

nominal 

2.074 

1.450 

1.450 

FLL 

2.099 

1.544 

1.544 


TABLE II: LDA+DMFT results for SrVOs at T = 200 K and 
U = 10 eV. Both t2g and eg orbitals are treated by DMFT. 



FIG. 1: (Color online) LDA+DMFT total density of states 
for SrV0 3 using three different DC potentials. Experimental 
photoemission is reproduced from Ref. l22l (parameters listed 
in table [TT]). 

n° = 1. The actual value of the DC potential Vd c differs 
for less than 0.1 U, which leads to almost indistinguish¬ 
able spectra on the real axis, and from the previously 
published results [9], hence we do not reproduce them 
here. We found a general trend in all materials stud¬ 
ied that the exact DC is somewhat smaller then given 
by FLL formula. For Ce, the Hartree contribution to 
DC potential is Vh = n/U ~ 0.99717, the semi-local 
exchange contribution is 14 ~ —0.485 U and LDA corre¬ 
lation is V c « —0.08817, hence the total DC potential is 
Vh+V x +V c « 0.42417, which is slightly smaller than FLL 
formula t/(n/ — 1/2) — J/2(nf — 1) ss 0.53317 or nominal 
formula 17(n° — 1/2) — J/2(n° — 1) = 0.517. It is interest¬ 
ing to note that the semi-local exchange used in LDA is 
quite different from the exact exchange value. The latter 
is only \Vf\ = Un/14: ss 0.071 U, a substantially smaller 
value then the semi-local exchange |I4| ~ 0.48517. This 
shows why DC within LDA+DMFT is so different from 
the Hartree-Fock value of the DMFT self-energy, i.e., 
£(w = oo). 

Next we present tests for SrV 03 , which is a metallic 
transition metal oxide with nominally single electron in 
the t2g shell. Near the Fermi level £4, there are mostly 
t2g states. The majority of eg states are above Ep, how¬ 
ever, due to strong hybridization with oxygen some part 
of eg orbitals also gets filled. There are two ways the 
DMFT method can be used here. In the first case, one 
can treat only the t2g shell within DMFT. The vast ma¬ 
jority of DMFT calculations for SrV03 were done in this 
way. In this case, all three DC potentials again give very 
similar results and the spectra is almost indistinguishable 
from previously published results in Ref. [TO} One can also 
treat dynamically with DMFT the entire d shell. This 


TABLE III: LDA+DMFT results for LaV0 3 at T = 200 K 
and U = 10 eV. Only t2g orbitals are treated by DMFT. 


LaV03(t2g+eg) 

nt.2g+eg 

nt 2 g 

n eg 

yalg/jj 

vi 9 /u 

vi 9 

exact 

2.444 

2.048 

0.397 

1.596 

1.599 

1.665 

nominal 

2.344 

2.032 

0.312 

1.458 

1.458 

1.458 

FLL 

2.706 

2.167 

0.540 

2.114 

2.114 

2.114 


TABLE IV: LDA+DMFT results for LaVOs at T = 200 K 
and U = 10 eV. Both t2g and eg orbitals are treated by 
DMFT. 


case is presented in Table [TT] and spectra in Fig. [Tj One 
can notice that the exact and the nominal DC give very 
similar while the FLL formula gives 14% larger nd- 
This is because the value of the DC potential is substan¬ 
tially larger (« 40 %) when using FLL as compared to ex¬ 
act case. It is nevertheless comforting to see that 40% er¬ 
ror in double-counting still does not leads to major failure 
of LDA+DMFT. We plot the spectra in Fig. [T] to show 
how this change in Vd c leads to shift of oxygen-p spec¬ 
tra relative to vanadium-d states. For the exact DC, the 
oxygen peak positions match very well with the experi¬ 
mentally measured spectra. The nominal valence is quite 
close to the exact spectra, while FLL formula leads to an 
upward shift of oxygen for roughly 0.6 eV, which is still 
relatively small compared to the difference in the double¬ 
counting potentials, which is Vf c LL — V// xact « 5.37 eV. 



FIG. 2: (Color online) LDA+DMFT total density of states 
for LaV0 3 using the three different DC formulas, (a) only 
t2g orbitals are treated by DMFT (b) both t2g and eg or¬ 
bitals are treated dynamically. Experimental photoemission 
is reproduced from Ref. (23) 

Next we present results for the Mott insulating ox¬ 
ide LaVC> 3 , which is solved in two ways, i) treating only 
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the t2g orbitals dynamically with DMFT, presented in 
Table III and Fig. [2^,, and ii) treating both t2g and eg 
with DMFT. In the first case, the valences are similar in 
all three double-counting formulas. The t2g occupancy is 
very close to nominal value 2. The exact double-counting 
is again smaller than given by FLL or nominal formula, 
which leads to a slightly larger splitting between oxygen- 
p and V-d states, i.e., slight upward shift of oxygen states 


in Fig. 
in table 


ig.[2k 1 
iblejl V[ 


In case ii) displayed in Fig. §> and tabulated 
where both the t2g and eg orbitals are treated 
by DMFT, the FLL formula dramatically fails, as it over¬ 
estimates the valence, i.e., n^ LL — n exact ~ 0.26. While 
the Mott gap does not entirely collapse, it is severely 
underestimated by FLL formula. The nominal valence, 
however, gives very similar results as the exact DC. This 
improvement of nominal DC as compared to FLL was 
pointed our in Refs. TO] and was found to hold not 
just in transition metal oxides but also in actinides 121 1 . 
The t2g occupancy rit 2 g in the nominal and exact DC is 
very close to nominal value of 2, equal to the scheme i) 
presented above. It is therefore not surprising that the 
spectra in Fig. Sb and Fig. are similar, with slight 
improvement compared to experiment when eg orbitals 
are also treated by DMFT. 

In summary, we presented continuum representation 
of the Dynamical Mean Field Theory, which allowed us 
to derive an exact double-counting between Dynamical 
Mean Field Theory and Density Functional Theory. The 
implementation of exact double-counting for solids shows 
the improved agreement with experiment as compared to 
standard FLL formula. Previously introduced nominal 
DC formula mm is in very good agreement with exact 
double-counting derived here. 
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